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ABSTRACT 

We derive the projected surface mass distribution Em for spherically symmetric mass 
distributions having an arbitrary rotation curve. For a galaxy with a flat rotation 
curve and an ISM disk having a constant Toomre stability parameter, Q, the ISM 
surface mass density E ff as well as Sm both fall off as Rr 1 . We use published data on 
a sample of 20 well studied galaxies to show that ISM disks do maintain a constant Q 
over radii usually encompassing more than 50% of the H I mass. The power law slope 
in T, g covers a range of exponents and is well correlated with the slope in the epicyclic 
frequency. This implies that the ISM disk is responding to the potential, and hence 
that secular evolution is important for setting the structure of ISM disks. We show 
that the gas to total mass ratio should be anti-correlated with the maximum rotational 
velocity, and that the sample falls on the expected relationship. A very steep fall off 
in E s is required at the outermost radii to keep the mass and angular momentum 
content finite for typical rotation curve shapes, and is observed. The observation that 
H I traces dark matter over a significant range of radii in galaxies is thus due to 
the disks stabilising themselves in a normal dark matter dominated potential. This 
explanation is consistent with the cold dark matter paradigm. 

Key words: galaxies: structure - galaxies: evolution - galaxies: spiral - galaxies: 
irregular - dark matter. 
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1 INTRODUCTION 

H i has long bee n the best tracer of Dark Matter (DM ) 
in galaxies (e.g. iBosmal Il98ll : Ivan der Hulst et all Il995l ). 
This is because it typically extends much further than 
the optically bright portion of a galaxy, in a fairly reg- 
ular disk. But this ability to trace DM seems uncanny: 
typically the projected DM surface densit y scales very 
well wi t h the measured H I surface densi t y jBosmal 198 it 
Sancisi 1983; Carignan & Bcaulieu 1989; Cari gnan et alj 
199d: ICarignan fc Puchel Il990allbl: I Jobin 



Broeilsl Il992l; iMeurer et al.1 119961 : iHoekstra et all ^ 
Mestell (|l963l ) had already shown that the flat rotation 
curve typically seen in the outer parts of spiral galaxies 
requires a surface mass density S(i?) oc fall-off in 

the disk, if that is where the dominant mass is located. 
This is not what is observed in the optical but close to the 
behaviour of H I disks. H I not only is a good dynamical 
tracer of DM but its distribution scales linearly with the 
DM in the outskirts of galaxies. This has prompted some re- 
searchers to posit that DM may be gaseous, perhaps in a disk 



confi gu ration (jP fennige r et al. | 1 1994 IPfc nnigcr & Co mbes! 
1994 iGerhard fc Silkl Il996l: iPfenniger fc Revazl 120051 ; 
Bournaud et al.ll2007l : iHessman fc ZiebartJ 2011 ). The scal- 
ing between DM and H I is also well explained by the Modi- 
fied Newtonian Dynamics (MOND) hypothesis in which the 
gravitational force law is modified; in a MOND analysis, 
H I i s the only significant mass in the outskirts of galaxies 
(e.g. ISanders fc McGaughl [20021) . Either explanation poses 
a problem for the standard Cold Dark Matter (CDM) sce- 
nario. A gaseous form of baryonic DM would be dissipative, 
whereas in the CDM scenario DM only interacts via gravity 
and so is non-dissipative. The MOND scenario requires no 
DM. 

We propose an explanation for the linear scaling be- 
tween H I and DM that is consistent with the CDM sce- 
nario. The ISM distribution in disks is configured to main- 
tain a uniform minimum stability over as much of the disk 
as possible. For a flat rotation curve, this will result in a sur- 
face density profile having the same form as the dominant 
mass. In §[2] we give the basis of our model. In §[3] we gather 
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recent H I rotation curve (RC) data to test this hypothesis. 
§ [4] discusses our results and presents our conclusions. 



2 THE STRUCTURE OF DYNAMICALLY 
STABLE GAS DOMINATED DISKS 

The stability of a purely stellar or purely gaseo us disk is 
a well stud ied prob l em st arting with the work of ISafronovl 
ljl960h and lToomrd l|l964h . For a disk to be stable against 
axi-symmetric perturbations, the support from a combina- 
tion of random motions and centrifugal forces must be larger 
than the gravitational attraction. This is expressed as a ra- 
tio, the "Toomre Q" parameter: 



Q 



tvGT, 



(1) 



where E is the mass density in the disk, a is the velocity 
dispersion, and k is the epicyclic frequency given by 



1 + 



RdV_\ 
V dRJ 



(2) 



with V being the rotational velocity at radius R. For a 
purely gaseous disk to be stable Q > 1, while a purely 
stellar disk requires Q > 1.07 = 3.36/-7T for stability. 
Unstable disks result in the formation of bars and spi- 
ral arms which gather the ISM, enhancing star formation 
efficiency and thus reheating the disk t hrough feedback 
llHohll 1 19711 : ISellwood fc Carlberd 1 1984 iDebattista etafl 
120061 1. The stability of a multi-component di s k, i.e. gas and 



stars, is more complex ( Jog fc Solomon|[l984 lRafikovll200ll ; 
iRomeo fc Wiegerg|2011 ). Here we are concerned primarily 



with the outer disk, where the ISM mass dominates. In those 
cases the single fluid Q given by eq. [1] is sufficient for our 
purposes. In some of the galaxies we analyse the stellar disk 
is important, and we show that the results are usually not 
significantly different over the radii we are concerned with 
when a multi-component disk analysis is employed. 

As done commonly, we consider a galaxy where the 
dominant mass has a spherical distribution of mass with 
density pm(R), which also contains an embedded "light" 
(low mass) gas disk having surface mass density E s . We 
take the projected total surface mass density Em(-R) to be 



E M (#) = 2Rp M (R). 



(3) 



This is the projection of a spherical shell onto a ring of 
the same radius at the equator. This is the correct defini- 
tion when we are concerned with the rotation curve (here- 
after RC, or V) of the dominant mass. Under our assump- 
tions, disk particles in circular orbits only feel the poten- 
tial of the mass interior to their orbit. We assume this ge- 
ometry precisely be c ause it allows some simple derivations. 
lOlling fc Merrifieldl |200Cf ) note that theory indicates that 
DM halos should be somewhat flattened spheroids, but that 
the limited observations do not clearly state what the flat- 
tening typically is. If the DM is in a disk, as might be ex- 
pected for gaseous form of DM, then the spherical approx- 
imation does not hold. As can be discerned from Fig. 1 of 
iHessman fc Ziebartl (|201lh . assuming a spherical geometry 
will cause us to underestimate Em(-R) at large R in this case. 
However, for typically assumed RC shapes and beyond a few 
core radii, the ratio of true to inferred Em(-R) changes slowly 



and only by ~ 10%. Hence the simplifying assumption of a 
spherical halo is not critical to our results. 

Assuming standard Newtonian gravity and pure circu- 
lar orbits, it is straight-forward to derive Em(-R) as a func- 
tion of the RC 



Em = 



d(RV 2 



2-kGR dR 
which becomes 

1 V /„„ dV 



Em = 



+ v 



2itG R 

We define £ to be the local gas to total mass ratio 
C| E„ 



(4) 



(5) 



(6) 



The case of £ being constant corresponds to the gas disc 
tracing the total mass. We will also have occasion to consider 
the integrated quantity 



Z 



M a 



M(R n 



(7) 



where M(i? max ) is the mass within the maximum radius, 
ij max , of the H I distribution. 

We emphasise that the mass densities Em, and pu in 
eq. [3] and throughout are the projected and spherically av- 
eraged mass densities of all mass, i.e. disk (stars and gas), 
bulge, and halo (luminous and dark). Likewise £ and Z are 
the local and integrated gas to total mass ratio. We adopt 
these definitions to simplify the analysis, avoiding the need 
to fit the rotation curve into contributions from the differ- 
ent components. The quantit y C is thus th e reciprocal of 
the ratio plotted in Fig. 7 of iBosmal (|l981T ) where it first 
became apparent that H I traces whatever dominates the 
mass in the outskirts of galaxies. The implications of our 
results stem from the well established result that at large 
radius the dominant form of ma s s in galaxies is DM (e.g . 
Freemarj|l970l; iRubin et al.lll978l: iFaber fc Gallaghe j Il979l: 



Kentl 19871: ICarignan fc Beaulieull 19891 : ICarignan et alJl99C 



de Blok et afll20Qg| ). 

Solving eq.[T]for E 9 and employing eq.[2]and eq.[5]yields 



c 



2a y/2[l + (R/V){dV/dR)] 



? V + 2R(dV/dR) 
It is instructive to adopt a power law form to the RC 
V w fcTT, 



(8) 



(9) 



where k is a constant. This approximation works fairly well 
over a limited range of radii in galaxies, which is sufficient 
for our purposes. This results in 



c = 



a V /8 ( 1 +7) 



(10) 



QkRt{2y + l)' 
The power law index 7 has a narrow range of allowed values: 
1 



< 7 < 1- 



(11) 



A 7 below this range means the RC drops faster than Ke- 
plerian, while 7 above this range requires a "hollow" mass 
distribution. 7 = —1/2 corresponds to £ = 00 since this re- 
quires Em = 0. In the outer disks of many disk galaxies the 
RC is flat at its maximum value, hence 7 = 0, V = Knax, 
and thus 
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V8 a 

Vmax Q 



(12) 



If a and Q are also constant, or their ratio is, then £ is also 
constant and the gas disk will track the total mass distribu- 
tion. From eq.[5]then 



1 V 

1 v ma; 

2tvG~ R~ 



E 9 (i?) 



tyQG~ W 



(13) 



This is the well known relationship of the total surface mass 
density falling off as _R _1 where the RC is flat, which is a fair 
approximation of what is observed in most spiral galaxies. 

We posit that disks evolve towards maintaining a con- 
stant Q. Simple feedback should encourage such a condition. 
Over the optical face of a galaxy, star formation is likely 
to be the regulating agent. Regions of the disk where Q is 
higher than average have a disk that is a combination of hot 
or under-dense compared to their surroundings. In these re- 
gions any star formation activity would decrease (relative 
to their surroundings), lowering a, thus allowing more ISM 
to accumulate or cool and decreasing Q. In regions where 
Q is low, the disk is a combination of cold or over-dense. 
Star formation will be enhanced in these areas, increasing a 
as feedback from the newly born stars kicks in. While the 
outer disk is usually considered to be devoi d of star forma- 
tion, the discovery of o utlying H H regions dFerg uson et al.l 
1 1998T I and XUV disks (jThilker et alj 120051 , 120071 ) indicates 
that in many cases there are sources of new stars that can 
help regulate disks. Even in the case o f pure gaseous disks, 
the simulations by IWada et all (2002) show that a turbu- 
lent clumpy disk develops with a large range of Q over short 
scales, but with a quasi-steady equilibrium maintained with 
little variation in average Q with time or radius. 



3 OUTER DISK Q MEASUREMENTS 

We gathered data on 20 galaxies to test our hypothesis that 
outer disks maintain a nearly constant Q. The majority of 
the data is from The H I Nearby Galaxy Survey ( THINGS) 
for wh ich detailed published RC s can be found in loh" et al.l 
(|2008h and lde Blok et all l|200St ). Data from individual stud- 
ies of six additional galaxies with extended disks are in- 



2009; Elson ct al. 



2010l ; IWarren et al 



eluded to test the robus t ness of the results llCannon et al 
poid: iGentile et all 120071: IStruve et al 
]20o3 ; IWestmeier et al.ll201ll v" The mair 



criteria for selection is that the galaxy have tabulated data 
available from recent studies (within ~ 5 years) and H I pro- 
files extending beyond the optical radius R25 (where the B 
band surface brightness is 25 mag arcsec -2 ). Table[T]lists ba- 
sic quantities of the sample and the data sources, arranged 
by the maximum rotational velocity Vmax- This sample cov- 
ers 50kms _1 <; Knax <, 375 km s" 1 and the full range of 
late-type galaxy morphologies including spirals from types 
SO (NGC 1167)' to Sd (e.g. NGC 300), irregulars of types Sm 
(IC 2574) and Im (e.g. NGC3741) as well as Blue Compact 
Dwarf (BCD) galaxies (e.g. NGC 2915). The majority of 
the sample is nearby, with only two having distance D > 15 
Mpc. Therefore, we adopt D values that are not based on 
the Hubble flow, where possible. The sources of D are given 
in Table [1] For the two cases where we use D based on 



redshift, we adopt the model given by NEE0 for the Hub- 
ble flow corrected for inflow to the Virgo cluster. Grea t At- 
tractor, and Shapley super-cluster (jMould et al.ll200oT ) and 
standard cosmological parameters Ho = 73 km s" 1 Mpc -1 , 

^matter 

= 0.27, and Q. 

vacuum — 0.73 

The main observational quantities of importance are 
the RC and the E g profile. The RCs are shown in Fig. [T] 
We fit a cubic-spline to the RCs using knots set by eye. 
In performing the fit, we keep the number of spline knots 
to a minimum and try to follow the data to within the er- 
rors. However, we smooth over small scale fluctuations in the 
RCs, presumably due to spiral arms or non circular motions. 
The fitted splines are shown as continuos lines in Fig. [1] We 
chose this functional form for its flexibility and because it al- 
lows for easy evaluation of the derivative dV/dR needed for 
the calculation of k and thus Q. For compariso n, we plot a 
model RC using the functional form adopted bv lLerov et al.l 
(2008), with parameters Va a t = 150 km s 1 , Ra at = 0.5i?25, 
and R25 = 8 kpc. The sample includes many galaxies with 
flat RCs at large radii, like this model (DD0154, ES02lS 
NGC 2915, NGC 2403, NGC 3198, NGC 6946, NGC 2841), 
or are still rising (IC 2574, NGC 3741, NGC 3621). How- 
ever, about half the sample have RCs that have a substan- 
tial range of radii where they are declining (NGC 2366, 
NGC 300, ADBS J113^, NGC 7793, NGC 4736, NGC 5055, 
NGC 2903, NGC 3521, NGC 3031, NGC 1167). To illustrate 
the range of shapes, we fit a power-law RC to the data be- 
tween a limited range of radii, defined below. These and 
the other power-law fits in this paper were performed as 
linear fits in log-l og space using the IDL procedure MPFIT 
l|Markwardtl 120091 ). The fits are shown as dashed lines in 
Fig.[T] with the fit parameters, including zeropoint, slope 7, 
and dispersion of the residuals e listed in Table 2. 

The E 9 profiles are shown in Fig. [2] They were derived 
from the inclination corrected H I profiles assuming E s = 
1.3Ehi to account for elements heavier than hydrogen. Since 
we are primarily interested in the outer disk we assume that 
the molecular content is negligible. The profiles are drawn 
in log-log space, to highlight any power-law behaviour. In 
general, the profiles are nearly flat or decreasing towards 
small radii, and steeply dropping at the largest radii, leaving 
a "knee" where the profile has an approximately power law 
form. We examined each profile and determined the inner 
and outer radii at which the power-law portion begins and 
ends, respectively, and the fraction of the total H I mass 
Mm within those radii. On average 21(±10)% of Mm is 
interior to the inner radius and 75(±10)% is interior to the 
outer, hence the power law knee contains a bit over 50% of 
Mm on average. In order to treat the galaxies consistently, 
we determine the radii enclosing 25% and 75% of Mm, Ri 
and R2 respectively, and fit a power law 

E 9 oc R N (14) 

to the points between these radii. The fits are shown in Fig. [2] 
with the parameters listed in the corners of the panels and 



1 The NASA/IPAC Extragalactic Database (NED) is operated 
by the Jet Propulsion Laboratory, California Institute of Tech- 
nology, under contract with the National Aeronautics and Space 
Administration. 

2 full name: ESO215-G?09 dWarren et alj|2004h. 

3 full name: ADBS J 113845+2008 llCannon et al.ll2009l) . 
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Figure 1. Rotation curves (RC), ordered and scaled by maximum rotational velocity, V m ax, plotted against radius R scaled by the 
maximum radius of the data -R max - We show the measured points and their errors (blue), the adopted spline RC fits (magenta), with 
the middle 50% of the HI mass highlighted (green), and a power law fit to the data in this region (dashed green). Each panel lists the 
galaxy name is at upper left, the optical radius R25 in kpc at lower left, V^n ax in km s — 1 at upper right and the power law slope 7 at 
lower right. Vertical (grey) lines mark intervals of R25, w hile horizontal lines mark intervals of 50 km s — 1 . The toy model plotted in the 
lower right panel is an RC of the form used by iLerov et al. (2008) with i?, flat = 4 kpc. 
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Figure 2. Radial profiles of ISM surface mass density E 9 (blue) plotted on a log-log scale. Line colour and highlighting are the same as 
in Fig.[T] Power law fits to the highlighted portion of the profile are shown with dashed lines (dark green). Each panel lists dispersion of 
the residuals of the fits at upper right, the power law index N at lower left, and log(S 9 ) at R25 at lower right. The toy model shown in 
the lower right panel is an exponential profile (eg. 1 1 5 1) with scale length a -1 = R25 and central density S s ,o = 10 M.q pc -2 . 
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Table 1. Sample properties 



Galaxy 
(1) 



morphology 



D 



(2) (3) 



R25 
(4) 



Mm 

(5) (6) 



Data ; distance source 



(7) (8) 



DD0154 

NGC3741 

ES0215 

NGC2366 

IC2574 

NGC2915 

NGC300 

ADBSJ1138 

NGC7793 

NGC2403 

NGC3198 

NGC3621 

NGC4736 

NGC6946 

NGC5055 

NGC2903 

NGC3521 

NGC3031 

NGC2841 

NGC1167 



IB(s)m 
Im 

Im a 
IB(s)m 
SAB(s)m 
BCD 6 
SA(s)d 
BCD C 
SA(s)d 
SAB(s)cd 
SB(rs)c 
SA(s)d 
(R)SA(r)ab 
SAB(rs)cd 
SA(rs)bc 
SA(rs)bc 
SAB(rs)bc 
SA(s)ab 
SA(r)b: 
SAO" 



1.3 
3.2 
5.2 
3.4 
3.8 
4.1 
2.0 

50.0 
3.9 
3.2 

13.8 
6.6 
4.7 
6.8 
7.9 
8.9 

10.7 
3.6 

14.1 

66.0 



1.24 
0.97 
0.80 
2.20 
7.26 
1.23 
5.81 
2.27 
5.78 
7.52 
13.22 
9.57 
5.41 
10.79 
14.30 
15.50 
11.49 
11.42 
28.90 
23.49 



8.28 

6.97 
10.82 

8.20 
11.07 
10.12 
19.36 
24.20 

7.74 
18.01 
37.51 
25.77 

9.61 
22.08 
38.15 
29.34 
31.17 
14.80 
51.68 
63.88 



0.45 
0.22 
1.46 
0.81 
1.67 
0.49 
2.55 
1.92 
1.11 
3.59 

14.50 
9.43 
0.52 
7.00 
8.29 
6.45 

12.70 
4.05 

13.90 

12.20 



50 
52 
52 
58 
78 
86 
99 
106 
118 
144 
159 
159 
198 
224 
212 
215 
233 
260 
324 
377 



de Blok et al. (2008) ; Makarova et al. (19981 
Gentile et al. (2007) ; Dalcanton et al. (2009) 
Warren et al. (20041 ; Karachentsev et al. (20071 
Dalcanton ct al. 



Oh et al. (2008 

Oh et al. (2008 1; Dalcanton et al 
Elson et al. (2010 1; Meurer et al. 
Wcstmeiei^e^iL^ (2011 1 ; Freedman et al 
Cannon et al 




dc Bk)k_et_al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
dc Blok ct al 
de Blok et al. (20081 
de Blok et al. (20081 
Struve et al. (20101 ; NED 




NED 

Karachcntsev^tji^ (2003b) 
Freedman et al. (2001 ) 
Freedman et al. (2001) 
Freedman et al. (2001) 

(2003a) 
(2000) 



Karachentsev ct al. 



Karachentsev ct al. 
Tullv et al. (2009) 
Drozdovsky & Karachentsev 

Freedman et al. (2001 ) 
Saha et al. (2006) 



Columns: (1) ga lax y name; (2) ga l axy m orpholog y from the RC3 catalogue (dc Vaucouleurs et al.lll99ll ) except as following: a. 
IWarren et al. I ll2004h . b.| Meurer et al. I lll994l ). c.[c annon et al. I d2009h ; (3) Distance in Mpc; (4) radius at the B band 25 mag arcsec 
(5) radius of the outermost H I measurement; (6) H I mass in 10 9 Mq', (7) Maximum rotation amplitude in km s~ 1 . 



in Table 2. The range Ri to R2 is highlighted in our figures 
by using thicker lines to plot the profiles. N spans the range 
-2.01 (NGC 300) to -0.32 (NGC 3621). The power law fit 
is reasonable. The average es g = 0.042 dex (10%), and in 
only two cases is e Es > 0.1 dex (in NGC 4736, NGC 5055). 

While we use a power law approximation of E 9 for con- 
venience, H 1 pr ofiles are not always char acterised as such 
in the literature. Ferguson fc Clarke! (|200ll ). for example, use 
an exponential profile: 

E 9 = E 9 , oe - afl (15) 

where E 9i o is the central ISM surface mass density and a -1 
is the disk scale length. They show that the ISM in their 
viscous disk models are well characterised by such a profile 
over a wide range in radii, with the ISM scale length greater 
than that of the stellar disk (which is well known to follow 
an exponential profile ; lFreeman|[l97bl ). The bottom right 
panel of Fig. [5] shows an exponential profile having E 9 ,o = 
10 Mq pc~ 2 and a" 1 — R25. This profile is similar in shape 
to that of our sample. We fit a power law to this profile 
in the same manner as the sample galaxies, and report the 
results in Table 2 and Fig [2] Over the relevant Ri to R2 a 
power law fit to an exponential profile gives a es of 0.03 
dex (7%). In comparison, 55% of our sample have es g less 
than or equal to this. Thus a power law is at least as good a 
functional form as an exponential for about half the sample. 

Figure[3]shows the £ profiles in order to test the premise 
that H 1 traces the total mass. To form the ratio we take E 9 
from the observed profile and T,m from eq. [5] In general the 
profiles are highly structured. This is due to irregularities in 
both the E 9 and V profiles. The ( profiles are particularly 
sensitive to the latter since it depends on the derivative of 
the RC (eq. [5]). The toy model, which combines the smooth 
RC shown in Fig. [T] and the smooth E 9 profile shown in 



Fig. [21 has = 0.07 dex. Six of the galaxies in the sample 
have £ profiles that have at this level or smaller. 

Figure U shows the Q profiles, for an assumed constant 
velocity dispersion a = 8kms _1 . The average log(Q) be- 
tween Ri and R2 and the e in this value are listed in Ta- 
ble 2 and the panels of Fig. 2] For the majority of the sam- 
ple Q is quite flat between R\ and R2, and somewhat be- 
yond for many cases. Over this range, the average eq is 0.06 
dex (15%), while in two cases it is more than twice that: 
NGC 4736 and NGC 5055, for which the variation in Q 
corresponds largely to irregular structure in the E 9 profile 
(Fig. [21 . For NGC 2366, Q becomes undefined for much of 
the disk beyond R2 because, in this case, its RC is declining 
more steeply than a Keplerian decline. Since this is unphys- 
ical, it is likely that a w arp or non-circular motions cause 
the RC o f loh" et alj (|2008l ) to be underestimated. In the bot- 
tom right panel we combine the toy models shown in Fig. [1] 
and Fig. [2] This model has eq = 0.09 dex, while 70% of our 
sample have values less than or equal to this. As expected, 
Q is rising beyond R2 for the majority of the sample. 

The a = 8kms _1 we adopt is in the middle of the pack 
compared t o wh at has been adopted in other studies; e.g. 
Kcnnicutt jl989l ) adopts a = 6 km s 1 , while iLerov et al.l 
(2008) use a = llkms -1 . Changing to a different constant 
a will on ly change Q by a con stant multiplicative value. 
Following iTamburro et"atl {2009) we also performed calcu- 
lations where a declined linearly with radius, using their 
profiles for the few overlapping cases between our study 
and theirs, and otherwise setting a = lOkms -1 at R25 and 
falling linearly with radius to a = 5 kins -1 at the last mea- 
sured point of the radial profiles. The resultant Q profile do 
not look very different from fig. [4] in particular, the decline 
in a does not remove the rise in Q often found at large R. 

While Fig. [4] shows that Q is fairly constant over radii 
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Figure 3. Profiles of the ratio £ of gas surface mass density E 3 to projected total surface mass density Em, plotted on a log scale, 
against linear scaled radius. Line colour and highlighting are the same as in Fig. [T] The average log(f) and dispersion about the average 
in the highlighted region is marked with the solid and dashed horizontal (gray) lines and listed at the lower left and lower right corners. 
The vertical lines show intervals of i?25- 
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Figure 4. Toomre stability parameter Q profiles on a log scale, plotted against linear scaled radius. Line colour and highlighting are the 
same as in Fig.[T] The average log(Q) in the highlighted region is marked at lower left and the rms about this average is given at lower 
right, and these are shown with the horizontal (grey) lines. The vertical lines show intervals of R.25- In the six cases where the fit range 
is largely within R25 we show profiles of the two fluid stability parameter Q21 (Romeo & Wicgcrt 2011) as red dotted lines (thicker over 
the fit range); the pink thick line shows the average log(Q2f) over the fit range and the pink dashed lines are offset from this by the rms 
about this average. 
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incorporating about half of the H I, usually it increases out- 
side this range. The explanation of what is happening at 
large radii is fairly simple. There, RCs typically are flat or 
become flatter as can be seen in Fig [T] As shown by eq. 1131 
for (" to be constant in this limit requires E s oc Such a 

profile can not be maintained to arbitrary radius because it 
would require infinite Mm- To have a disk with finite mass 
and angular momentum requires a more rapid drop-off in 
E 9 , as is observed beyon d R2 (Fig. El), hence Q typically 
rises at the largest radii. IHoekstra etall |200ll) point out 
that the sharp drop off in E 9 at large R does not corre- 
spond to an expected decline in RCs, thus casting doubt on 



the ability of H I to trace DM. iHessman fc Zicbart 



counter that the fitting technique of IHoekstra et al 



2011 



2001 



over emphasises the RC fits at large R where much of the 
hydrogen is ionised. Neutral or ionised, there is a finite ISM 
mass, and this will limit the ability of the ISM to trace DM. 

There are multiple causes for the rising Q profiles to- 
wards small radii. The central parts of galaxies typically 
have a significant molecular content, which we have not in- 
cluded here, hence Q is overestimated. In addition, the stel- 
lar disk usually dominates the mass distribution, so that a 
single compon ent Q is inadequate for determining the tru e 
disk stability l|Lerov et al. 1 120081 ; [r omco fc Wiegertll20Tl . 
Our assumption has been that the H I largey resides in the 
outer disk. However, in six of the sample galaxies (IC 2574, 
NGC 7793, 4736, 2903, 3031, and 2841) the fitting range is 
largely interior to R25 , hence the stellar disk is likely to play 
an important role in determining the stability of the disk in 
these cases. 

In order to determine the true disk stability for these 
cases, we calculated the two fluid (stars and gas) stability pa- 
rameter Q2f- There are various formulations of the two fluid 
stabi l ity parameter fe.g.lJog fc Solomonll 19841 ; IWang fc Silkl 
1 19941 ; iRafikovllioOll ; iRomeo fc WiegertlbOlIfT Here, we cal- 
culate Q2i from eg. 9 of I Romeo fc Wiegertl (|201ll ) which 
accounts for the thickness of the stellar and gaseous disks. 
It can be rewritten in a simplified form as 



1 

Q2I 



+ 



Pa 



(16) 



Where and Q a are the stability parameters for the stel- 
lar and gaseous components of the disk calculated separately 
using eq. [TJ and p*, p g are weight factors which depend on 
the velocity dispersions of the stars and gas and the value 
of Qi, compared to Q g . Calculation of Q21 requires the stel- 
lar mass density E* in the disk which we derive from the 
stellar mass radial profiles of lde Blok et al.l (|2008l ). and the 
radial component of the stellar velocit y dispersion el l ipsoid 
o\t,r- We take cr*, r ~ cr* jZ /0.6 following IShapiro et alj (|2003l ) 
where <t* jZ is the vertical component of the stellar velocity 
dispersion. This is given by 



o"*, z = \J 2-7rGEdisk/i* 



(17) 



where Edisk = E* + E 9 is the total disk surface mass density 
and h i, is the stellar disk scale height. Following iLerov et al.l 
(|200St ) we estimate the scale height from the disk scale length 
using hi, ~ Z*/7.3 and fit the exponential scale length U over 
the same fit region as highlighted in our figures. While we are 
primarily are concerned with Q^i over the fitted region it is 
instructive to also see its behaviour beyond this range. In the 
central regions, in cases where the measured E* is greater 



than the extrapolated exponential fit, we replace E* with 
this fit, under the assumption that the excess light represents 
a bulge or thick disk with a scale height larger than the disk 
and thus has a lesser contribution to the disk potential than 
is expected for the surface brightness. Similarly, at large R, 
beyond the last measured E* , we also adopt the extrapolated 
fit when calculating Q2f ■ The results of the Q2i calculations 
are shown in Fig [4] as the red dotted lines. The average 
log(Q2f) over the fitted range is shown as the thick pink 
lines and the dashed pink lines are offset from this line by 
its rms. 

Figure [4] demonstrates that the Q2t profiles are flatter 
than the Q = Q g profiles when the entire radial range is 
considered. The biggest difference is for NGC 4736 where 
the average ± dispersion is (log(Q2f)) — 0.53 ±0.06 over the 
fitted region; the dispersion is a factor of three lower than 
that for the gas Q. For the remaining cases, the difference be- 
tween Q and are Q2t more subtle, over the fitted range where 
we find <log(Q 2f )> = 0.17 ± 0.05, 0.46 ± 0.02, 0.052 ± 0.04, 
0.30 ± 0.04, and 0.59 ± 0.12 for IC 2574, NGC 7793, 2903, 
3031, and 2841 respectively. The dispersion in log(Q2f) is 
smaller than that of log(Q) for half the cases: NGC 4736, 
2903, and 3031. In most cases Q is an adequate proxy for 
the total disk stability parameter Q2f at the radii we are 
interested in here, even when much of the H I is within the 
optical radius. However, a more sophisticated analysis, such 
as using Q2t is required to extend the analysis even closer 
to the center, or in cases like NGC 4736 where the stel- 
lar disk strongly d ominates the gaseous disk at all radii. In 
IZheng et all (|2012t ) we perform such an analysis over the op- 
tically bright portion of galaxies, using various prescriptions 
for Q2t- There we show that the assumption of a constant 
stability disk can explain the relative distributions of gas, 
stars and star formation over much of the optically bright 
portion of galaxies. 

Close examination of Fig. [T] and Fig. [2] shows that an- 
other of our premises is not exactly correct: the RCs are 
not always flat. Instead there is a rather wide distribution 
of rotation curve power law indices 7. What then causes the 
nearly constant Q seen in Fig. This requires the numera- 
tor and denominator in the defining equation (eq.[T} to have 
the same shape, i.e. the same slope. For a constant cr, this 
requires the k and E 9 profiles to have the same shape. To 
test this, we fit a power law k oc R m to the epicyclic fre- 
quency profile. The k profiles and the fits to them are shown 
in Fig. [S] with the relevant fit parameters listed in Table 2 
and the panels of Fig. [5] A power law form is a reasonable 
approximation to the k profiles between _Ri and i?2- 

Figure [6] plots the power law index in E 9 , N against 
the power law index in k, M. There is a crude correla- 
tion between the indices; the Pearson's correlation coeffi- 
cient r xy = 0.45, with a 2% chance that the correlation is 
random. Examination of the figure shows that there are two 
obvious outliers, and we can see plausible reasons for their 
discrepant behaviour in each case. NGC 3031 (M81) is in a 
nearby highly interactive group with three close companions 
(M82, NGC 3077, and HoIX). These may affect the outer 
ISM distribution of NGC 3031 through either stripping ma- 
terial, or having material stripped from them. NGC 300 has 
a very extended and lopsided H I distribution and the steep- 
est N — —2.01 in our sample. IWestmeier et all l|201ll ) note 
that there are morphological signs of ram-pressure stripping 
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Figure 5. Epicyclic frequency k profiles on a log-log scale. Line colour and highlighting are the same as in Fig.[T] In each panel, the rms 
of the power-law fits (green dashed lines) are listed at upper right, the power law index M at lower left, and log(K) at H25 at lower right. 
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Figure 6. Index N of the fit to the E 9 profiles in Fig. [2] plotted 
against the index M of the fit to the re profiles shown in Fig. \E\ 
The diagonal line marks N = M. Outliers are plotted as hollow 
circles, are labelled, and discussed further in the text. 

of the outer disk, which could steepen the E 9 profile. They 
also note that their ATCA data may not capture the to- 
tal H I flux due to missing short spacings. The E 9 profile 
we use was de rived from a data cub e that combined the 
ATCA data of IWestmeier et all i|201ll ) and single dish H I 
data fro m the Parkes 64m tel escope obtained for the GASS 
project (jKalberla et al-l fioiO), thus recovering the H I flux 
missing from the ATCA observation^. Excluding these two 
galaxies r xy — 0.63, with a 0.3% chance of being random. 
We conclude that there is modest correlation between the 
indices which scatter about the line N = M, that is the k 
and E s profiles have the same shape. This is exactly what 
is needed for a constant Q disk. 

A further test of the hypothesis that galaxy disks evolve 
to constant Q is given by eq. [T^] which predicts that £ and 
Knax should be anti-correlated. These quantities are plot- 
ted (in the log) in Fig. [7] The filled circles show the aver- 
age Iog(C) (i.e. (log(£))), between Ri and R2 as marked in 
fig. [3] The expected anti-correlation is clearly present, with 
r xy = —0.60; the probability that this is due to random sam- 
pling of uncorrelated data is less than 0.3%. The solid line 
shows the expected anti-correlation for galaxies with a flat 
RC, eq. I121 and for the assumed a = 8kms _1 and Q = 3 
(this corresponds to the average log(Q) for our sample). The 
dashed line is a % 2 fit to the data: 

(log(C)) = (0.83 ± 0.13) - (0.93 ± 0.06) log(V max ). (18) 

The dispersion about this fit ei og (c) = 0-28 dex, while the 
average offset between the expected relation and this fit is 
0.08 dex, very close to what is expected from the standard 
error on the mean, 0.06 dex. This demonstrates that eq. [12] 
provides a reasonable model for the average gas to total mass 
ratio between Ri and R-z. 

The hollow circles in Fig. [7| show the ratio Z plotted 

4 Without the Parkes data we would have an even steeper TV = 
-2.36. 
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Figure 7. Average log(f) between Ri and R2, (log(f)), plotted 
as filled circles against rotation curve maximum Vmax. The er- 
ror bars represent e^. The hollow circles represent the integrated 
ratio Z for the same galaxies. The solid line shows the expected 
relationship for galaxies with a fiat RC, eq. I12I and a = 8kms _1 , 
Q = 3. The dashed line shows the fit to the (log(f)) data re- 
ported in eq. I18I while the dotted line shows the fit to the Z data 
reported in eq. I19I 

against V max . These quantities are even more strongly anti- 
correlated with r xy — —0.80 and a less than 0.002% chance 
that the correlation is due to random sampling of uncorre- 
lated data. The dotted line shows an equally weighted least 
squares fit to the data 

\og(Z) = (1.17 ± 0.19) - (1.12 ± 0.09) log(V max ); (19) 

the dispersion about this fit £iog(c) — 0.24 dex. This is 
tighter than eq. 1181 probably because Z is better defined 
than {log(f)) which is more susceptible to noise in the E 9 
and V profiles. 



4 DISCUSSION AND CONCLUSIONS 

The correlation between the shape of the k and E 9 profiles 
shown in Fig. [6] is profound. Since M is determined purely 
by the potential, this implies that the ISM disk is responding 
to the potential, and hence that secular evolution is driving 
the correlation. An alternate hypothesis is that the structure 
of galactic disks is set by the m a ss and angular mom entum 
accretion history l|Barnesl [20021 ; i ancisi et alj I2008T ). How- 
ever, it is not clear why there should be any such corre- 
lation under this scenario. CDM simulations indicate that 
Milky- Way mass galaxies h ave had typically only two major 
merg er event since z ~ 2 l|D'onghia et al.ll2008l ; I Cole et aU 
l2007h. 10 Gyr ago, which agrees with observations of merg- 
ers (fconselice et al.ll2008T h So, to the extent the accretion is 
from major-mergers, the outer disks should be very lumpy. 
Our sample does include some lumpy disks (e.g. NGC 4736 
and NGC 5055) and also some galaxies that may currently 
be interacting (NGC 300, NGC 3031) perhaps contributing 
to the scatter and outliers in Fig. [6] The fact that the major- 
ity of the sample falls on the N = M line suggests that most 
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are not strongly interacting and that prior interactions hap- 
pened far enough in the past that the disk has re-stabilised 
and smoothed out to trace k. The timescale for doing this 
is the orbital time i or b (~ 3tdyn where idyn is the dynami- 
cal timescale). For our sample, the average t OT h{R2) = 0.6 
Gyr. Taking the separation between major merger events 
to be 3 Gyr, there should be approximately five orbits be- 
tween mergers, sufficient time for structure in the disk to 
smoot h out. Alternatively, if accretion is s low "cold accre- 
tion" (|Keres et al.ll2005l ; ISancisi et al"1l2008l ) the disk would 
not be expected to be lumpy. However, that scenario does 
not provide an obvious explanation for the E g and k pro- 
files following each other. Of course, cold accretion combined 
with feedback to equalize Q, as advocated here, is consistent 
with our results. 

Previous studies have no ted the large scatter in the 
HI to total or DM ratio (e.g. lBosma|[l98ll ; iHoekstra et al.l 
l200lf ) . This is hard to explain in the context of H I be ing a 
linear tracer of a baryonic DM l|Hoekstra et al.ll200ll ). We 
show that an inverse correlation between £ and Vinax is ex- 
pected from constant Q disks (eq. [12) , and indeed is ob- 
served (Fig. [7). This avoids the need to place implausibly 
large quantities of molecular gas in disks. By linking a lu- 
minous component of the disk to the rotational velocity the 
Z - Knax relation is r eminiscent of the Tully-Fisher relation 
iTully fc Fished Il977j) and the baryonic T ullv-Fisher rela- 
tion (|Freemanlll999l ; lMcGaugh et alj|2000h . However, since 
Z is a ratio, it does not provide a means to estimate ISM 
mass from V mxl . To do so, one needs to know the extent 
of the H 1 mass distribution, since in the case of 7 = the 
local £ and the integrated gas to total mass ratio within R 
remain constant and equal to each other. Alternatively, if 
one has Mm and Kiax one can estimate the extent of the 
ISM distribution. This is typically the case in blind H I sur- 
yeys, such as HIPASS l|Barnes et all 120011 ), and ALFALFA 
(|Giovanelli et al.ll2005h where the H I flux is known but the 
source is unresolved. Then one can use eq. [19] to estimate Z 
and from that the maximum H I extent 



Rn 



GMi 



(20) 



This may be useful for determining whether follow-up ob- 
servations of a particular galaxy are likely to be fruitful, 
or estimating the covering factor of H I absorbers (e.g. 
IZwaan et ai1l2005l ). 

Our results imply that secular processes are i mpor - 
tant for setting disk structure. iLvnden-Bell fe Pringlel (ll974T ) 
noted that viscous disks should evolve so that mass is con- 
centrated in the centre and angular-momentum goes to infin- 
ity. As dissipative encounters cause ISM mass to be funnelled 
towards the galaxy centre, the disk must also spread to 
conserve angular mo mentum, and the ISM di sk size should 
grow with time (e.g. iFerguson fc~C larkc 200l]). High resolu- 
tion simulations of galactic disks show that transient den- 
sity waves can increase the size of disks and alter their 
metallicity distributions (Roska r et al. 2008a,b). The flat- 
tening of metallicity gradients l|Werk et al.l |2010| . 1201 lh in 
outer disks may be evidence of ISM circulation in spread- 
ing disks. If disks are spreading now, they should have been 
denser and more compact in the past. This would imply 
a higher molecular fraction from the increased hydrostatic 
pressure, as well as increased total gas content (to account 



for the stars that have since formed). iBraun et al.l (|201lh 
noted an order of magnitude increase in the molecular mass 
for the most luminous star forming galaxies from now to 
z ~ 0.4, mirroring the incr ease in the cosmic star for mation 
rate density psfr(z) (e.g. iHopkins" fc Beacomll2006l ). They 
argue that this implies that recent star formation evolution 
is largely due to th e run down in the available ISM supply. 
lHanish et all (120061 ) made a similar argument based on the 
slope in Psfr(- z ) being similar to that expected from the star 
formation law. As outer disks spread and evolve, Q should 
remain constant with radius. The exact values of Q and a 
will be set by feedback (with less efficient star formation in 
the thinning disk) and the available angular momentum. 

Finally, our results do not disprove MOND, nor do they 
rule out the possibility that some or most of the DM is in a 
gaseous form. However, some of the appeal of these the- 
ories is that it was not clear under the prevailing CDM 
paradigm why dissipative gaseous disks should trace the 
non-dissipative CDM halo which resides in a spheroid. Our 
work provides this linkage by showing that ISM disks will 
trace DM as a natural consequence of disk stabilisation, and 
our tests of the H I dominated outer disks are consistent 
with that interpretation. 
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Note added in proof: We note a previous paper that 
reached similar conclusions that had escaped our attention. 
Struck-Marcell (1991, ApJ, 368, 348) used a related ap- 
proach, the assumption that gas disks maintain a balance 
of hyrdostatic forces, to derive a disk structure of E 9 oc R~ x 
for disks having a flat RC. That paper briefly alludes to the 
relationship with DM, but does not generalize the problem 
to arbitrary RC shape as we have done. 
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Table 2. Fitted quantities 
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Columns: (1) gives the galaxy name; (2) and (3) give the radii in kpc enclosing 
25% and 75% of the H I mass respectively; columns (4)- (6) relate to the power 
law fits to the rotation curves - (4) gives the zeropoint, (5) the power-law index, 
and (6) the dispersion of the residuals in dex; columns (7)- (9) relate to the power 

■ law fits to the E 9 (ISM surface mass density) profiles, with the individual columns 
following the pattern of (4)-(6); columns (10)-(12) relate to the power law fits 

■ to the k (epicyclic frequency) profiles, with the individual columns following the 
pattern of (4)-(6); columns (13) and (14) give the average logarithm of £ (local 
gas fraction) between R\ and R2 and its dispersion, respectively; columns (15) 
and (16) give the average logarithm of Q (disk stability parameter) between 
Ri and R2 and its dispersion, respectively. Values in parenthesis are the formal 
uncertainties in the power law slopes from the fits to the relevant profiles. 



